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We present a general method which expresses a unitary operator by the product of operators 
allowed by the Hamiltonian of spin-1/2 systems. In this method, the generator of an operator 
is found first, and then the generator is expanded by the base operators of the product operator 
formalism. Finally, the base operators disallowed by the Hamiltonian, including more than two- 
body interaction operators, are replaced by allowed ones by the axes transformation and coupling 
order reduction technique. This method directly provides pulse sequences for the nuclear magnetic 
resonance quantum computer, and can be generally applied to other systems. 
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I. INTRODUCTION 

In 1973, Bennett [|] proposed a reversible Turing ma- 
chine that is as efficient as an irreversible one, and this 
led to the idea of using a quantum system as a computer 
because the time evolution of a quantum system is re- 
versible. Feynman Q introduced the concept of a quan- 
tum computer, and its theoretical model was given by 
Deutsch [|. On the other hand, Fredkin and ToffoH |] 
proved that an arbitrary computation can be performed 
by a reversible Turing machine by showing that and, 
OR, and NOT gates can be generated by reversible 3-bit 
gates among which a Toffoli gate [|| is most frequently 
used nowadays. In quantum computation, a three-bit 
gate cannot be implemented directly because it requires 
a simultaneous interaction of three particles. Thus there 
have been efforts to find two-bit universal gates • In 

particular, Barenco et al. showed that a combination of 
two-bit c-NOT gates and one-bit gates can replace a Tof- 
foli gate, and proposed a method to make general n-bit 
controlled gates Therefore, it is proved that an arbi- 
trary computation can be performed by a quantum com- 
puter, and the implementation of these universal gates 
became the basic requirement for any quantum system 
to be a quantum computer. 

However, the proof that an arbitrary computation can 
be done by a quantum computer does not necessarily 
mean that we know a general implementation procedure. 
If n unitary operator U, equivalent to a combination of 
gates, is related to the Hamiltonian 7i of a certain quan- 
tum system hy U ~ cxp{—tHt/h), it can be reahzed by 
the time evolution of the system during time t. But there 
are only a few operations that can be implemented in this 
way by the limited Hamiltonians of nature. Therefore, it 
is very necessary to find a general method to implement 
an arbitrary operation using only the given Hamiltoni- 
ans. Feynman proposed a way to construct an artifi- 
cial Hamiltonian when U is given hy U — Uk ■ ■ ■ U3U2U1 
and all Hi's corresponding to Ui's exist in nature [Q, 
but it is impractical to construct artificial Hamiltonians. 
It will be more practical to partially control a Hamilto- 
nian by turning perturbations "on" and "off" if U can 



be expressed as a product of operators corresponding to 
the perturbation terms. Whether Feynman's artificial 
Hamiltonians or switchable perturbations are used, an 
operator of interest must be expressed as a product of 
the operators allowed by Hamiltonians. This is equiva- 
lent to finding the combination of universal gates, or a 
quantum network, and generally is a very difficult prob- 
lem having several solutions. 

In this work, we propose a general method of express- 
ing a unitary operator as a product of the operators al- 
lowed by the Hamiltonian of the spin-i systems including 
the nuclear magnetic resonance (NMR) quantum com- 
puter. This method makes use of the fact that a unitary 
operator U is always given hy U = exp[— iG] , where G is a 
Hermitian operator. Once the generator of an operator, 
G, is found, it is expanded by suitable base operators. 
Then U is expressed as a product of operators having 
only one base operator as a generator and, finally, each 
operator in the product is replaced by the allowed ones. 
Compared to a previous report , only the operators of 
physical variables were used in each transformation pro- 
cedure. This helps to understand the physical meaning 
of operations done by a quantum computer. 

II. DECOMPOSITION 

The first step of implementation is to find the genera- 
tor of a given operator. Since the only way to implement 
an operator is to use the time evolution of a state under a 
suitable Hamiltonian, a generator, which is the product 
of Hamiltonian and time, gives the physical information 
necessary for implementation. A unitary operator is rep- 
resented by a normal matrix and always diagonalized by 
unitary transformation. The matrix T that diagonalizes 
U also diagonalizes G as 

f/' = rC/Tt = e-'^«^' =e-^<=', (1) 

where U' and G' are diagonalized matrices of U and G, 
respectively. Once the operator and its generator become 
diagonal, G' is easily obtained from 
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-^G' 



(2) 

and G is obtained by inverse transformation G = T'^G'T. 
Since G is Hermitian, G"^^., the eigenvalues of G, are real 
and U'^.^. are complex with absolute value of unity. It is 
worthwhile to note that the mapping from C/^^, to G'^.^, is 
not unique. 

To relate the generator G with Hamiltonians, consider 
the following operators of the product operator formalism 
for N spin-i particles 



(3) 



where s = {ai, a2, a^} and is 0, x, y, or z. Iq is E, 
i.e., a 2 X 2 unity matrix, /q,. is a spin angular momentum 
operator for ^ 0, and q is the number of nonzero a^'s. 
For example, {Bs} for = 2 is given by 



q = 0: E/2 



hx , Ily , hz, hx , hy , hz 

'^hxhx, '2hxhyi '^hxhz, ■■■ 



(4) 



which are 16 Dirac matrices except the factor of i. In 
Eq. (^) , unity matrices are not shown and spin indices are 
added for convenience. {Bg}, consisting of 4^ elements, 
makes a complete set and, therefore, an arbitrary 2^ x 2^ 
matrix can be expanded by the linear combination of 
Bg^s. Since G and S^'s are Hermitian, coefficients of 
the linear expansion are real numbers and obtained by 
applying the inner product of G and i?s's. 

A unitary operator is now expressed as [/ = 
exp{—iJ2sbsBs), of which the generator is related to 
physical observables. In general, there exists no Hamil- 
tonian that corresponds to a linear combination of i?s's. 
Therefore, our next step is to express U as a product 
of single operators, which have only one Bg as a gener- 
ator like exp[— i6s_Bs]. Sometimes this decomposition is 
the most difficult step, and it has not yet been proven 
whether the decomposition is generally possible even for 
spin operators. Fortunately, many useful gates can be 
easily decomposed by using the commutation relations 
of -Bs's. Bs's are either commuting or anticommuting 
with each other. If G is expanded with only commuting 
-Bs's, U can be easily represented by a product of single 
operators as 



U = cxp[— z bgBs] ]^exp[— 



(5) 



A swap gate and an /-controlled phase shift gate used in 
Grover's search algorithm belong to this case. 

Even though a generator has non-commuting i?s's, 
there are cases where decomposition is straightforward. 
Suppose two base operators B^i and Bs2 satisfy the re- 
lation {h = 1): 

[BsuBs2]=iB,s] (6) 

then Bs3 also belongs to {Bg}. This commutation rela- 
tion makes the three operators B^i, Bs2, and Bsj, trans- 
form like Cartesian coordinates under rotation, meaning 
that 



exp[-i(/>i?s3]i?si(exp[-i0Ss3])''' = Bsi cos(/> -f Bs2 sin(/> 

(7) 

for cyclic permutations of si, s2, and s3. If a genera- 
tor has only these operators, it can be decomposed using 
Euler rotations. For example, ey.p[—i(j){Bsi +i?s2)] is un- 
derstood to be the rotation with the angle of about 
the axis 45° off the "i?si axis" on the plane of Bgi and 
Bs2 axes. Therefore, this operation is equivalent to the 
successive rotations about B^i and Bg^ axes as follows: 



-i0(S,i+Bs2) 



(8) 



This decomposition technique by Euler rotations is also 
applicable when an operator has a generator in the fac- 
torized form as follows: 



U — exp 



N 



n E< 



(9) 



where (jjia^ are real numbers. Since hx^ hy, and hz sat- 
isfy the commutation relation in Eq. and commute 
with any other spin operators with i ^ 1, spin-1 compo- 
nents are decomposed as 



Ui exp 



N 



yio 



E- 



^iho.,)Y[ 



i=2 



(10) 



where Ui is the product of the single operators of which 
the generators have only spin-1 components, correspond- 
ing to Euler rotations. Repeated applications of this pro- 
cess to successive spins give 



U = UN---Uie-'^Ul---U 



t 

AT' 



where 



N 



G = l[i<PME + cj),h^J. 



(11) 



(12) 



Then decomposition is finished because all terms in 
Eq. (|l^) commute with each other. All the controlled 
gates belong to this case. If none of the above methods 
are applicable, U can be approximately expanded as a 
product of single operators to any desired accuracy 19|. 



III. REDUCTION 

Although Bs is a product of spin operators that are 
physical quantities, not all Bs's exist in Hamiltonians 
that nature allows. The next step of implementation is 
to replace disallowed single operators in the product by 
allowed ones. The Hamiltonian of a spin-i system used 
for implementation of a quantum computer allows only 
the following single operators in general. 
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Ria{4>) 
Jija{4>) 



-l<plia, 

-i<p21ialj. 



(13) 



The first term is a rotation operator that rotates spin i 
about the a axis by the angle of 0, and the second one 
is a spin-spin interaction operator between spins i and 
j. The angle 4> in the second term is proportional to the 
spin-spin coupling constant and evolution time, but we 
denote it as a rotation angle because the effect of spin- 
spin interaction can be understood as a rotation of one 
spin due to the magnetic field of the other. Before going 
further, we assume the following more restricted set of 
operators as allowed ones in this study: 



(a 



= X OT y), 



(14) 



In this set, only x and y axes are used for single-spin ro- 
tations and a spin-spin interaction is limited to the Ising 
type. Needless to say, the greater the number of single 
operators allowed, the easier it is to implement an al- 
gorithm. However, Eq. (|l^ ) is a sufficient set to realize 
any unitary operators as shown below, and in fact these 
are the only operators allowed by an NMR quantum com- 
puter. Single-spin rotations are implemented by selective 
rf pulses and spin-spin interactions by Hamiltonian evo- 
lution with intermediate refocusing pulses |20| . These 
two rotation operators can generate any single bit opera- 
tion and the interaction operator can make a c-not gate 
in combination with rotation operators |^,|l|. There- 
fore, these three operators consist the minimum set to 
implement universal gates. 

Now, we are to show that the minimum set in Eq. ( |l^ ) 
can generate all the other single operators. First, the- 
single bit operation excluded in Eq. (p^), Ri.^{(j)), can be 
transformed from Rix{4>) as 

RU^) = R^y (-|) R^A^)R^V (|) • (15) 

This is the composite pulse technique well-known in the 
NMR experiments ||l^ . Any rotation about one axis can 
be replaced by the composite of rotations about the other 
two. All the second-order operators, where the nth-order 
operator means the single operator that has a generator 
Bg with q = n, can be transformed into the Ising-type op- 
erator in Eq. ( [T^ ) by this technique |^,^. For example, 
exp[— 102/^2; /j^] is transformed as 



ex.p[-uf>21ixlj^] = R, 



ly 



(16) 

In the same way, any nth-order operator can be trans- 
formed into the product of single operators and the nth- 
order coupling operator that is defined as the nth-order 
operator with all aj = z. 

After all the spin coordinates are changed to z using 
this technique, the operators with more than two-body 



interaction can be reduced to an Ising-type two-body in- 
teraction operator as discussed below. The key idea of 
the coupling order reduction is that the nth-order cou- 
pling operator can be thought as the (n — l)th-order one 
controlled by one spin state. For example, the third-order 
coupling operator [Q, e-xp^-icjj'ilizl j yJkz] is represented 
by 



exp[-«04/j^/j^/fe2] 



exp 







exp[-Z(/)2(/^ (g)/^)] 


Jjfe(0) 

Jjk{-i 






exp[i(/)2(/2 (g) 



(17) 



in the subspace of spin i. The final form of Eq. ( p_7D im- 
plies that the third-order coupling operator can be un- 
derstood as a second-order one with coupling between 
spin j and k, but its rotation direction depends on the 
state of spin i. We note that if one spin is flipped dur- 
ing the evolution of the spin-spin interaction, then the 
sign of the interaction changes and this has the effect 
of time reversal. This means that the rotation direction 
changes [ ^7|j2C| ] and, therefore, we can implement Eq. ([l^ ) 
with the second-order coupling operator by flipping spin 
i or k depending on the state of spin i. It is a well-known 
c-NOT (xor) gate that flips one spin depending on the 
state of the other spin. A c-not gate is given by 



NOT — 



R.. (I) R,. (I) Rjy (f ) J^, {-\) R,y (-f 



(18) 



up to an overall phase, and this is the product of allowed 
operators in Eq. (p^. 

In the same way, an nth-order coupling operator can 
be reduced to an (n — l)th-order one by conditionally 
flipping odd number of spins except the spin i. Repeated 
applications of this process obviously reduce an nth-order 
coupling operator to a second-order one. Fig. ^ shows 
the quantum networks of the nth-order coupling opera- 
tor and its equivalent combination of the allowed opera- 
tors. In Fig. |l](b), the c-not gates after the second-order 
coupling operator are inserted to flip spins to their orig- 
inal states. Instead of the c-not gates before and after 
the second-order coupling operator, a pseudo c-not gate 
Uij — Rjx{^)Jij{^)Rjy{^) and Ujj can be used, respec- 
tively. 

As an example, we apply this general implementation 
procedure to a Toffoli gate 1 25 1 . The generator of a Tof- 
foli gate obtained after the processes of diagonalization 
and inverse unitary transformation is expanded by base 
operators as 



j21izlsx 



Since all terms in this generator commute with each 
other, the corresponding operator is easily expressed by 
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the product of single operators. After replacing disal- 
lowed operators, hzhx and hzhzhx in this case, by 
allowed ones by axes transformation and coupling order 
reduction, the gate is finally expressed as 

Rlz{j)R2z{j)Jl2{-j)R3x{j)R3v{^)Jl3{-j) 
XJ23(-f)C/l2J23(f);7|2^3,(-f) ^ ^ 

up to an overall phase. 



IV. CONCLUSION 



Our method is applicable to any quantum computers 
that use spin-i states as qubits, because the operators 



of Eq. ( |14| ) make the minimum set required to those 
quantum computers. This method can be generalized for 
other quantum computer systems that should provide a 
complete set of operators similar to those of Eq. (14). 
Since our method uses generators that are closely related 
to Hamiltonians, this method helps us to see the phys- 
ical meaning of an operation. Operators with genera- 
tors disallowed by Hamiltonians are replaced by allowed 
ones using axes transformation and order reduction tech- 
niques. Therefore, it is possible to simulate an Hamilto- 
nian that does not exist in nature using this method, 
including more than two-body interactions. 

This method does not necessarily give either optimal 
or unique solution to implementation. As the number of 
spins increases, the number of base operators grows ex- 
ponentially, and a generator could have too many terms. 
Therefore, it is impractical to apply this method to a 
system with many spins, but our method still provides a 
good guide to implement an operator of interest. 



[13: 



[14 



[is: 
[16: 



[17: 
[is: 

[19 
[20 
[21 
[22: 
[23: 

[24: 

[25: 



A. Barenco, C. H. Bennett, R. Cleve, D. P. DiVincenzo, 
N. Margolus, P. Shor, T. Sleator, J. Smolin, and H. We- 
infurter, Phys. Rev. A 52, 3457 (1995). 
R. P. Feynman, in Feynman Lectures on Computation, 
edited by A. J. G. Hey and R. W. Allen (Addison- Wesley, 
New York, 1996). 



R. R. Tucci, preprint, |quant-ph/9902062 . 

O. W. Sorensen, G. W. Eich, M. H. Levitt, G. Boden- 

hausen, and E. Ernst, Prog. NMR Spectrosc. 16, 163 

(1983). 

E. Ernst, G. Bodenhausen, and A. Wokaun, Pnnciples 
of Nuclear Magnetic Resonance in One and Two Dimen- 
sions (Oxford University Press, Oxford, 1987). 
S. S. Somaroo, D. G. Cory, and T. F. Havel, Phys. Lett. 
A 240, 1 (1998) 

A. T. Sornborger and E. D. Stewart, Phys. Rev. A 60, 
1956 (1999). 

N. Linden, H. Barjat, J. Carbojo, and R. Freeman, 
Chem. Phys. Lett. 305, 28 (1998). 

J. Jones, R. H. Hansen, and M. Mosca, J. Magn. Reson. 
135, 353 (1998). 

T. F. Ha vel, S. S. Somaroo, C .-H. Tseng, and D. G. Cory, 
preprint, ^uant-ph/98 12086 

M. D. Price, S. S. Somaroo, C. H. Tsendg, J. C. Gores, 
A. F. Fahmy, T. F. Havel, and D. G. Cory, J. Magn. 
Reson. 140, 371 (1999). 

Experimental simulation of the third order coupling op- 
erator was shown in a different way by C. H. Tseng, 
S. S. Somaroo, Y. Sharf, E. Knill, R. Laflamme, 
T. F. Havel, and D. G. Cory, Phys. Rev. A 61, 012302 
(1999). 

Recently, a different method to construct n-bit Toffoli 
gates was reported by M. D. Price, S. S. Somaroo, 
A. E. Dunlop, T. F. Havel, and D. G. Cory, Phys. Rev. 
A 60, 2777 (1999). 



[1] 
[2] 
[3] 
[4] 

[5] 



[6] 
[7] 

[8] 

[9] 
[10] 

[11] 
[12] 



C. H. Bennett, IBM J. Res. Dev. 6, 525 (1973). 

R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 

D. Deutsch, Proc. R. Soc. London Ser. A 400, 97 (1985). 

E. Fredkin and T. ToffoU, Int. J. Theor. Phys. 21 219 
(1982). 

T. Toffoli, in Automata, Languages and Programming, 
edited by J. W. de Bakker and J. van Leeuwen (Springer, 
New York, 1980). 

D. Deutsch, Proc. R. Soc. London Ser. A 425, 73 (1989). 
D. P. DiVincenzo, Phys. Rev. A 51, 1015 (1995). 
D. Deutsch, A. Barenco, and A. Ekert, Proc. R. Soc. 
London Ser. A 449, 669 (1995). 

A. Barenco, Proc. R. Soc. London Ser. A 449, 678 (1995). 
S. Lloyd, Phys. Rev. Lett. 75, 346 (1995). 
A. Barenco, D. Deutsch, A. Ekert, and R. Jozsa, Phys. 
Rev. Lett. 74, 4083 (1995). 

T. Sleator and H. Weinfurter, Phys. Rev. Lett. 74, 4087 
(1995). 



-e- 



-e- 



(a) 



(b) 



FIG. 1. Quantum network for the nth-order coupling op- 
erator (a) and its equivalent network consisting of allowed 
operators (b). 
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